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Abstract 


This study aims to assess the intraspecific genetic diversity of Ixora macrophylla, a widespread species of Philippine 
Ixora recorded from several islands of the Philippines, and Z. auriculata, an endemic species that has yet to be included in 
molecular studies of the Philippine Zxora. The number of haplotypes, haplotype diversity, and nucleotide diversity of 19 
trnL-F and 17 ITS sequences of Z. macrophylla, as well as 4 trnL-F and 4 ITS sequences of I. auriculata were obtained 
using DNAsp 5.10.1 software, and pairwise distances were calculated using the nucleotide Kimura 2-parameter using 
MEGA 6.06. Two haplotypes of trnL-F and nine haplotypes of ITS were identified in six populations of Z. macrophylla. 
One of the two haplotypes of trnL-F was unique to the Mindoro population. One of the nine haplotypes of ITS was 
common among seven individuals from four populations. Intraspecific pairwise distances ranged from 0 to 0.1% for trnL- 
F and 0 to 0.9% for ITS. Mantel test showed weak correlations between the genetic and geographic distances for both 
trnL-F (r = -0.0380) and ITS (r = 0.0980) sequences. For the genetic diversity of Z. auriculata, two haplotypes of trnL-F 
and four haplotypes of ITS were identified, with intraspecific pairwise distances ranging from 0 to 0.1% in trnL-F and 0.3 
to 3.5% in ITS. The results for genetic diversity may be used to better understand the population genetics of the Philippine 


Ixora and provide insights for conservation. 
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Introduction 


Genetic diversity is represented by the gene pool, or the 
variation in genes, within and among populations (Rao & 
Hodgkin, 2001), and is influenced by factors such as natural 
selection, mutations, gene flow, and genetic drift (Amos & 
Harwood, 1998). Investigating the genetic diversity of species 
can reveal the patterns of gene flow among existing populations 
(Mason-Gamer et al., 1995). Geographical range is also an 
important factor in high genetic diversity among populations, as 
geographical isolation restricts the occurrence of gene flow 
among populations (Maki et al., 2003). A high population 


‘Department of Biological Sciences, College of Science University of 


Santo Tomas, España Boulevard, 1015 Manila, Philippines 

?The Graduate School, University of Santo Tomas, España Boulevard, 
1015 Manila, Philippines 

Research Centre for the Natural and Applied Sciences, University of 
Santo Tomas, Espafia Boulevard, 1015 Manila, Philippines. 


*Corresponding email: cibanag@ust.edu.ph 


Date Submitted: 31 July 2019 
Date Accepted: 04 May 2020 


© Association of Systematic Biologists of the Philippines 


genetic diversity in abundant species within populations can be 
associated with larger population sizes, and high rates of gene 
flow with other populations (Hague & Routman, 2016). Genetic 
diversity has significant implications to conservation and 
speciation, as it is linked to the adaptability of species and their 
populations to cope with environmental stresses (Frankham, 
2005), and can drive the evolution of species through natural 
selection. With climate change being a primary concern for 
ecologists, the capability of species to adapt to environmental 
changes is a necessity for their survival. 

The Philippines is considered a leading biological hotspot, 
as it supports a significant number of endemic flora and fauna 
(Myers et al., 2000) and, as an archipelago, serves as a valuable 
model for studies concerning phylogenetics and population 
genetics (Brown et al., 2013). Among the diverse flora are the 
endemic species of the pantropical genus [xora L. With 
approximately 530 species worldwide and 30 out of 41 
Philippines species are endemic (Alejandro, 2007), Ixora is the 
third largest genus under the Rubiaceae which grows 
predominantly in low- to mid-altitude humid rainforests (Davis 
et al., 2009). With its tropical nature, the Philippines provides a 
suitable environment for the propagation of such species. Recent 
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Figure 1. Jxora macrophylla (Cagayan) -(A) stipule, (B) infructuscence, (C) cauliflorous inflorescence (D) 
inflorescence with open flowers. Ixora auriculata (Bulacan)- E. stipule, (B) infructuscence, (C) terminal inflorescence 
(D) inflorescence with open flowers. Photos taken by Sarah Grace S. Zamudio (A-D), Cecilia I. Banag (E-G), and Pieter 
B. Pelser & Julie F. Barcelona (H). Scale bars: 1mm (A-F,H); 10mm (G). 


Philippine Journal of Systematic Biology Online ISSN: 2508-0342 Volume 14 Issue 1 - 2020 | 2 


Casipit et al.: Genetic diversity of Ixora 


studies on the molecular phylogeny and biogeographical 
patterns of Philippine xora showed that its polyphyly is 
attributed to at least five colonization events, of which one 
lineage is derived from species from Wallacea, and the 
remaining four lineages related to species from Asia. The 
incongruence between the topologies of nuclear and chloroplast 
datasets characterized the complex evolutionary history of the 
Philippine Zxora. Furthermore, the close relationship between 
widespread species and locally adapted narrow endemics 
suggests that the widespread species serve as pools from which 
speciation takes place (Banag et al., 2017). 

Ixora macrophylla Bartl. (Figs. 1A-C) is endemic to the 
Philippines and is widely distributed in most islands and 
provinces of Luzon, Visayas, and Mindanao. This species can 
be distinguished from other Philippine /xora primarily through 
its cauliflorous inflorescences, a unique feature in the genus, 
whose inflorescence is usually terminal. Ixora auriculata Elmer 
(Figs. 1D-G) is endemic only to several provinces in Luzon. It 
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Figure 2. Map of the Philippines with marked localities where the 
utilized samples were collected. 


© Association of Systematic Biologists of the Philippines 


is characterized by its sessile leaves with cordate or auriculate 
bases, pendulous inflorescences, and white corollas. 

Since I. macrophylla have a wider distribution in the 
Philippines compared with other single island endemics, often 
their conservation is not a major concern for policy makers. 
However, these widespread species can provide opportunities to 
investigate the population genetics of endemic Jxora in the 
Philippine archipelago. The purpose of this research is to 
provide baseline information on the genetic diversity of Z. 
macrophylla among six populations namely; Cagayan, 
Pangasinan, Aurora, Marinduque, Occidental Mindoro, and 
Surigao del Norte. The chosen populations are representative 
throughout the Philippines to maximize geographical distances. 

In the recent molecular study on Philippine xora, some 
endemic species such as Z. auriculata were not included in the 
phylogenetic analysis due to lack of available leaf material. 
Recent field collections conducted in two localities; Angat, 
Bulacan and Sibuyan Island, Romblon by the Thomasian 
Angiosperm Barcoding and Phylogeny Group (TAPBG) 
resulted in the recollection of Ixora species identified as Z. 
auriculata based on morphology. However, after thorough 
examination of the reproductive parts of the samples from Angat 
and Sibuyan, several differences mostly in the color and 
measurements of the inflorescence axes were observed raising 
the question whether the samples from these localities are 
conspecific or not. To date, molecular data on I. auriculata have 
yet to be produced. 

Thus, the objective of this study is to assess the 
intraspecific genetic diversity of the six Z macrophylla 
populations and the two J. auriculata populations using the 
nuclear ribosomal marker internal transcribed spacer (ITS) and 
the chloroplast marker trnL-F. 


Materials and Methods 


Taxon sampling 

Fieldworks were conducted in different provinces of the 
Philippines (Fig. 2) after securing Gratuitous permits and Local 
transport permits from the Department of Environment and 
Natural Resources (DENR), City Environment and Natural 
Resource Office (CENRO) and Provincial Environment and 
Natural Resource Office (PENRO) to collect plant samples of J. 
macrophylla and I. auriculata. The fruits and flowers were 
preserved in 70% ethyl alcohol for storage and leaf samples 
used for molecular analyses were preserved in silica gel (Chase 
& Hills, 1991). The herbarium specimens were mounted onto 
herbarium sheets and deposited at the University of Santo 
Tomas Herbarium (USTH) for accessioning. Refer to Appendix 
1 (List of herbarium specimens examined). 
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Table 1. Primers used in the amplification of nuclear and chloroplast regions 


Region Primer Primer Sequence (5’ -> 3’) Reference 

ITS P16F TCACTGAACCTTATCATTTAGAGGA Popp & Oxelman, 2001 
P25R GGGTAGTCCCGCCTGACCTG 

trnL-F C CGAAATCGGTAGACGCTACG Taberlet et al., 1991 
F ATTTGAACTGGTGACACGAG 


DNA Extraction, Amplification, Purification, and Sequencing 

Leaves stored in silica gel were homogenized into powder 
using an Omni Bead Ruptor. Genomic DNA was extracted from 
dried leaf samples using the DNeasy Plant Mini Kit (Qiagen, 
Germany) by following the protocol recommended by the 
manufacturer. The desired markers were amplified through 
polymerase chain reaction (PCR) using the Biometra T Gradient 
thermocycler. Both the ITS and ¢trnL-F gene regions were 
amplified from extracted DNA of I. macrophylla and I. 
auriculata (Table 1). The amplicons were subjected to agarose 
gel electrophoresis to confirm the presence of DNA bands. The 
amplified PCR products were purified using QlIAquick 
Purification Kit (Qiagen, Germany) following the protocol 
provided by the manufacturer. The purified DNA samples were 
then sent to Macrogen Inc. (South Korea) for bidirectional 
sequencing. The sequences were assembled using CodonCode 
Aligner v7.1.1 (CodonCode Co. USA), and ambiguous bases 
were manually edited based on readings of the 
electropherogram. The sequences were manually aligned using 
Mesquite v3.51 (Maddison & Maddison, 2018). 


Statistical analyses 

For the genetic and statistical analyses, the ITS and trnL-F 
sequences of I. macrophylla generated in this study were added 
to the sequences of 7. macrophylla that were collected from 
Aurora, Occidental Mindoro, and Surigao del Norte used in the 
preliminary study conducted by Gaticales et al. (2016). The 
number of haplotypes (h), haplotype diversity, (Hd) and 
nucleotide diversity (x) within and among populations of 7. 
macrophylla was determined using DnaSP version 5.10.1 
(Librado & Rozas, 2009). Pairwise distances of 7. macrophylla 
were calculated using the nucleotide Kimura 2-parameter model 
using MEGA version 6.0 (Tamura et al., 2013). The same 
analyses were conducted with the two populations of J 
auriculata collected from Bulacan and Romblon. 

The correlation between genetic and geographic distances 
of I. macrophylla populations was evaluated using the Mantel 
test. For the first matrix, pairwise genetic distances of the 19 
individuals from the 782 base pairs (bp) of trnL-F and 17 
individuals from the 674 bp for ITS were computed as 
implemented in MEGA v6.06. The geographic distances 
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between the collection sites were computed with an Excel macro 
for the Vincenty formula (Dalgleish, 2015), obtaining a second 
matrix. Correlations and the Mantel test with 999 iterations were 
calculated with the ExtraStats function in Poptools version 3.2.5 
(Hood, 2010). Mantel test was not done for I. auriculata study 
since only two populations were included in this study. 


Results 


Genetic diversity of I. macrophylla 

This study used a total of 19 trnL-F sequences and 17 ITS 
sequences obtained from six different populations of Z. 
macrophylla. Of these, eight trnL-F and eight ITS sequences are 
newly generated in this study. Two haplotypes were identified in 
the trnL-F sequences, while nine haplotypes were found in the 
ITS sequences. One of the haplotypes for trnL-F was common 
among all the individuals of five populations: Cagayan, 
Pangasinan, Aurora, Marinduque, and Surigao, while the other 
haplotype was restricted to all samples from the Mindoro 
population. One of the haplotypes for ITS was common among 
seven individuals overall: two individuals from Cagayan, two 
from Marinduque, one from Mindoro, and two from Surigao. 
Another haplotype is common among all individuals of the 
Pangasinan population. The rest of the individuals, particularly 
the three from Aurora, one from Cagayan, two from Mindoro, 
and one from Surigao all possess their own unique haplotypes. 

The average haplotype diversity for the trnL-F sequences 
was 0.28070, while the average haplotype diversity for the ITS 
sequences was 0.82353. There is no haplotype diversity for the 
trnL-F sequences within populations. On the other hand, ITS 
sequences are generally more genetically diverse, ranging from 
0 (Pangasinan and Marinduque) to 1 (Aurora and 
Mindoro). The nucleotide diversity for the trnL-F sequences is 
0 within all populations yet yielded a value of 0.00036 among 
the populations. Nucleotide diversity of the ITS sequences 
ranged from 0 (Pangasinan and Marinduque) to 0.00405 
(Aurora), with an overall value of 0.00284. The data for the 
genetic diversity of the six Z. macrophylla populations are 
summarized in Table 2. The pairwise distance of the Z. 
macrophylla sequences ranged from 0 to 0.1% in trnL-F, and 0 
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to 0.9% in ITS. The smallest value of 0 was usually observed 
between samples within populations. 


Mantel Test 

The Mantel test between the genetic and geographic 
distances of Z. macrophylla resulted to a weak negative 
correlation coefficient of -0.0380 for trnL-F and a weak positive 
correlation coefficient of 0.0980 for ITS (Table 3). As 
represented by the diamonds in the charts (Fig. 3), both values 
for the correlation coefficient fall within the 95% confidence 
interval of -0.1668 to 0.2528 (trnL-F) and -0.3413 to 0.3516 
(ITS). 


Genetic diversity of I. auriculata 

A total of four trnL-F and four ITS sequences of Z. 
auriculata were generated. Two haplotypes were retrieved from 
the trnL-F sequences, while four haplotypes were identified 
from the ITS sequences. One of the two haplotypes of trnL-F 
was common among two individuals from Bulacan and one 
from Romblon, while the other haplotype was unique to one 
individual from Romblon. On the other hand, all four 
individuals possess unique haplotypes for ITS. The average 
values of haplotype diversity for both datasets are 0.5 for trnL-F 
and | for ITS. The haplotype diversity for the trnL-F sequences 
was 0 for the Bulacan samples, and 1 for Romblon samples, 
while the haplotype diversity for the ITS sequences yielded a 
value of 1 for both the Bulacan and Romblon samples. The data 
are summarized in Table 4. 

The pairwise differences of the Z. auriculata sequences 
ranged from 0 to 0.1% in trnL-F and 0.3% to 3.2% in ITS. The 
mean genetic distances of ITS were 0.45% within populations 
and 3.05% between the two populations. For trnL-F, the 
smallest value of 0 was observed among the two individuals of 
Bulacan and one individual of Romblon, while the value of 
0.1% was obtained between one of the Romblon individuals and 


Table 2. Genetic diversity of the six 7. macrophylla populations 
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Figure 3. Chart representation of the Mantel test results for trnL-F 
(left) and ITS (right) sequences of 7. macrophylla 


all other individuals. For ITS, smaller values (0.3-0.6%) were 
observed between individuals of the same population, while 
larger values (2.9-3.2%) were detected between individuals from 
different populations. 


Discussion 


Genetic diversity of I. macrophylla 

Haplotypes, or sets of alleles that are inherited together 
from a single parent, can be used in inferring patterns of gene 
flow among populations (Templeton, 2004). The number of 
haplotypes for the trnL-F region of I. macrophylla is relatively 
lower compared to that of their ITS. Among all 19 individuals, 
16 share a common haplotype, which include the samples from 
Cagayan, Pangasinan, Aurora, Marinduque, and Surigao. This is 
evidence that gene flow had occurred among these five 
populations. The other haplotype is restricted to the three 
samples from Mindoro. The low haplotype diversity of trnL-F 
can be associated with its less variation and lower substitution 
rates (Andreasen et al., 1999; Dong et al., 2012). Despite its 


A. trnL-F B. ITS 

Group BOE, maes diey diees, laity Noot Noon “Guerin, Nt 
Cagayan 3 1 0.00000 0.00000 0.00000 3 2 0.66667 2.00000 0.00304 
Pangasinan 3 1 0.00000 0.00000 0.00000 3 1 0.00000 0.00000 0.00000 
Aurora 4 1 0.00000 0.00000 0.00000 3 3 1.00000 2.66667 0.00405 
Marinduque 2 1 0.00000 0.00000 0.00000 2 1 0.00000 0.00000 0.00000 
Mindoro 3 1 0.00000 0.00000 0.00000 3 3 1.00000 2.00000 0.00304 
Surigao 4 1 0.00000 0.00000 0.00000 3 2 0.66667 0.66667 0.00101 
Total Data 19 2 0.28070 0.28070 0.00036 17 9 0.82353 1.86765 0.00284 


© Association of Systematic Biologists of the Philippines 


Volume 14 Issue 1 - 2020 | 5 


Casipit et al.: Genetic diversity of Ixora 


Table 3. Results of the Mantel test for 7. macrophylla populations 


A. trnL-F 

Correlation between 

matrices 

r= -0.0380 

Randomized Matrices 

(999 runs) 

rmin= -0.1857 

r min95 -0.1668 95% 
confidence 
interval, 
min 

r ave -0.0021 Average 

r max95 0.2528 95% 
confidence 
interval, 

r max= 0.5821 max 


naturally low substitution rate, a unique character in the trnL-F 
sequences is observed within the population of Mindoro. The 
restricted haplotype of the Mindoro population may be 
attributed to the geologic history of Mindoro, as the island was 
separated from the rest of the Philippines during the Pleistocene 
period (Heaney, 1986 as cited in Jones & Kennedy, 2008). The 
isolation of Mindoro from the rest of the populations studied 
during the Pleistocene period may have enabled the Mindoro 
population to evolve independently, which may have led to its 
restricted haplotype. 

The genetic diversity of 17 individuals across the six 
populations can be associated with the nine identified 
haplotypes in the ITS sequences. This is supported by the 
recombination and higher substitution rates that occur within 
nrDNA, thus increasing the variation among individuals, and 
allowing them to evolve quickly (Alvarez & Wendel, 2003). 
The populations of Pangasinan and Aurora possess their own 
unique haplotypes. This suggests the lack of apparent gene flow 
for these populations, since they do not share common 
haplotypes with other populations. High haplotype diversity and 
low nucleotide diversity suggest that the populations recently 
diverged from each other, before the time period long enough to 
observe mutations that resulted to different haplotypes, but not 
long enough to accumulate mutations for bigger differences 
among sequences (Garg & Mishra, 2017; Su et al., 2015). 

For the ITS sequences, representatives from Cagayan, 
Marinduque, Mindoro, and Surigao share the same haplotype, 
which indicates that gene flow had occurred among these 
populations. This is contrary with the results of the number of 
haplotypes for the trnL-F region of I. macrophylla, suggesting 
that Mindoro had not exchanged genes with other populations 
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A. ITS 

Correlation between 

matrices 

r= 0.0980 

Randomized Matrices 

(999 runs) 

rmin= -0.5345 

r min95 -0.3413 95% 
confidence 
interval, 
min 

r ave -0.0155 Average 

r max95 0.3516 95% 
confidence 
interval, 

r max= 0.6051 max 


and the haplotype acquired from gene flow is not shared among 
the individuals from Mindoro. This inconsistency may be 
attributed to the incongruence between chloroplast and nuclear 
datasets of the Philippine xora. Chloroplast DNA, which is 
maternally inherited, is used to explain the biogeography of the 
genus, as taxa with similar haplotypes may have originated from 
the same location. On the other hand, the biparentally-inherited 
nuclear DNA reflects the relatedness between taxa in an 
evolutionary sense (Banag et al., 2017). 

The Pangasinan and Marinduque populations have shown 
to possess the least genetic diversity. This makes them more 
susceptible to inbreeding depression, which lowers the overall 
fitness of the populations. In addition, populations with lower 
genetic diversity are less capable of responding to 
environmental stresses. Therefore, these populations should be 
prioritized for conservation. They can be rescued by inducing 
anthropogenic gene flow into these populations to increase their 
genetic variation and reduce the occurrence of inbreeding 
depression (Waller, 2015). 


Mantel test 

The Mantel test is a statistical test that determines the 
correlation of genetic distance with geographical distances. 
Mantel test showed that these two variables are weakly 
correlated. The correlation coefficient, as represented by the 
diamonds, falls within the range of statistical confidence 
intervals (Fig. 3). Therefore, the genetic distances for the six 
populations studied cannot be attributed to isolation-by-distance, 
in which populations that are geographically closer are 
genetically more similar than populations that are farther apart 
(Mims et al., 2016). This is expected due to gene flow between 
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Table 4. Genetic diversity of the two Z. auriculata populations 


A. trnL-F B. ITS 
haplotype Ave. no. of . R haplotype Ave. no. of ; 
Group No.of No.of diversity differences, Nucleo tide No. of No. of diversity differences, Núcleo tide. 
sequences haplotypes Hd K diversity, pi sequences haplotypes Hd K diversity, pi 
Bulacan 2 1 0.00000 0.00000 0.00000 2 2 1.00000 4.00000 0.00622 
Romblon 2 2 1.00000 1.00000 0.00127 2 2 1.00000 2.00000 0.00311 
Total Data 4 2 0.50000 0.50000 0.00063 4 4 1.00000 13.66667 0.02125 


populations, as inferred from the previous section (Feng et al., 
2014). The correlation between the two variables is expected to 
increase by adding more samples from different localities to this 
study. 

Aside from geographical range, other factors, such as 
dispersal capabilities, can influence the distribution of 
intraspecific genetic differentiation and gene flow. This is made 
possible by the long-distance seed dispersal by birds and 
mammals (Cain et al., 2000). Gravity seed-dispersal by birds 
and mammals has shown to be a factor in the genetic variation 
for Ammopiptanthus Cheng f. of the Leguminosae (Ge et al., 
2005). Drupes can promote long-distance dispersal, as they may 
be the most specialized fruit type for bird dispersal (Bremer & 
Eriksson, 1992). Frugivorous birds and small mammals that 
consume drupaceous fruits of Ixora can serve as dispersal 
vectors for the genus (De Block, 1998). 


Genetic diversity of I. auriculata 

The genetic diversity of the Z. auriculata samples is 
represented by two different haplotypes of trnL-F and four 
distinct haplotypes of ITS. In agreement with the previous 
section, the lower number of haplotypes identified in trnL-F 
relative to that of ITS sequences is attributed to their 
substitution rates. A common haplotype for trnL-F is shared 
between the two populations, thus inferring that the two 
populations may have originated from the same mother 
population (Banag et al., 2017). Four haplotypes of ITS were 
identified for Z. auriculata, with each haplotype unique to each 
individual. Therefore, the two populations of I. auriculata, at 
least for the samples utilized in the study, are genetically 
diverse. The absence of gene flow can be a major factor for the 
genetic differentiation between two populations (Feng et al., 
2014). 

The values for pairwise distances between the two 
populations of I. auriculata are notably high, with a mean 
distance of 3.05%, as compared with the pairwise distances 
observed among 7. macrophylla populations, whose distances 
peaked at only 0.9%. This raises the question whether the two 
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samples are conspecific. Several studies have utilized high 
genetic differentiation as one basis in determining whether 
samples are conspecific or not (Alejandro et al., 2007; 
Goodwillie & Stiller, 2001; Shneyer & Kotseruba, 2015). The 
high degree of genetic differentiation between the two 
populations of I. auriculata may suggest that these two are 
separate species, but the data presented are insufficient to 
support this claim. Phylogenetic analyses and more descriptive 
morphological examination of the two populations are necessary 
to reinforce this hypothesis. 


Conclusions and Recommendations 


The results of this study may serve as baseline information 
for further studies on Z. macrophylla, and I. auriculata and may 
be useful in designing appropriate conservation strategies. It is 
recommended that more samples from other localities are 
collected to validate the genetic diversity of these species, and to 
re-evaluate the correlation of the genetic and geographic 
distances of the Z. macrophylla samples. The high degree of 
genetic variation between the two populations of Z. auriculata 
may suggest that the two samples are separate species. Further 
studies should focus on performing phylogenetic analyses and 
providing detailed morphological descriptions and comparisons 
of the two populations to support this claim. 
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